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1. Introduction. 


Two active areas of research in spectral methods are the definition of well- posed 
approximations to the Navier-Stokes equations in which rigorous error bounds can be obtained, 
and the construction of spectral-type techniques which are applicable to problems in general 
domains. As regards the former, recent progress has established "staggered-mesh" formulations 
resulting in well-posed, solvable schemes that are optimal in the velocity [27], [4], [5], [22], 
[2]. As regard spectral techniques for simulation of flows in complex domains, several algorithms 
have been proposed, including multidomain collocation patching schemes [11], [23], [20], [21 ], 
[ 15] and spectral element variational techniques [25] , [13], [26], 

The first schemes which involved staggered meshes, resulted in well-posed problems (i.e. , 
without any spurious mode in the pressure) for one-dimensional problems [27], [4], however 
the extension of the ideas to higher space dimensions [22] introduced spurious modes into the 
system [2]. Although a workable scheme can be achieved by filtering the pressure [22], it is 
clearly desirable, in particular in three space dimensions, to construct a method in which the 
problem is intrinsically well-posed. To this end, a collocation technique is proposed in [2] , in 
which the velocity and pressure spaces are chosen so as to give a unique solution. For this last 
scheme, an error analysis has been performed, and spectral accuracy is proved. 

The spectral element spatial discretization involves a variational projection operator 
applied to elemental tensor product Lagrangian inter polants through local Chebyshev [25] , [13] 
or Gauss-Lobatto Legendre collocation points [26]. The technique is capable of handling general 
geometries with relative ease, due to the "automatic" patching inherent in the variational 
formulation. In this paper, we consider a synthesis of the staggered-mesh and spectral element 
concepts that represents an enhancement of both ideas. The spectral element discretization benefits 
by a significant improvement in the treatment of the pressure as compared to past methods, in 
which either spurious pressure modes are present [25], or in which the pressure is treated 
inaccurately [ 1 3], The staggered mesh schemes are improved in that the error estimates obtained 
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for the variational spectral element discretization are better than those obtained previously for 
collocation. 

In Section 2, we present our numerical method as applied to a Legendre spectral 
element-Fourier discretization of the steady Stokes problem. All results presented here extend 
directly to higher space dimension, as well as to the unsteady case, as will be discussed in future 
papers (see e.g. [ 1 9]). In Section 3 a theoretical analysis is performed, in which it is shown that 
no spurious modes appear in the pressure. Furthermore, optimal error estimates are obtained for 
both the velocity and pressure. Lastly, in Section 4, we present some numerical results. These are 
in accordance with the theoretical estimates, We provide also some details on the numerical 
implementation of the method that uses a new algorithm for solving spectral Stokes 
discretizations, the details of which will be presented in a future paper [ 1 8], 

In what follows, for any integer m and any domain A in DR or OR 2 , we denote by C m (A) the 
space of all functions that are continuous over A as well as all their derivatives up to the order m , 
and by C~( A) the space of all functions that are infinitely differentiable with compact support in 
A. 

In order to precise the sense in which the equations we shall consider have to be understood, 
we introduce some functional spaces. We denote by L 2 (A) the standard Lebesgue space provided 
with the norm || . || A and the scalar product (,,.) A ( or || , || and (...) when no confusion can occur), 
and for any positive real number r, the usual Sobolev space H r (A) provided with the norm || . || r A 
and semi-norm | . | r A (or || . || r and | . | r when no confusion can occur). 

Finally, C will stand for various constants that may vary from one line to the other, and for 
any function f depending on one variable x, we denote by f x the derivative of f with respect to x. 

Acknowledgments : The authors want to thank C. Bernardi for helpful comments concerning this 
work. 
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2. The Numerical Method, 


2J-T.he m od el problem 


Throughout the paper, Q is the domain Ax© , with A = ]-1 ,1 [ and © = ]0,2n[. The 
generic point in 0 will be denoted bg x = (x.y). We consider the Stokes problem in the domain 0 
for the velocity u = (u,v) and the pressure p, with no-slip boundary conditions in the first 
direction and periodic boundary conditions in the second. The problem is : Find (u,p) such that 


( 2 . 1 ) 


-vAu + Vp = f in 0 , 


[ div u = 0 in Q , 
with the following boundary condition 

(2.2) V y € © , u(- 1 ,y) = u( 1 ,y) = 0 , 

(2.3) V x € A , u(x,0) = u(x,2n) . 

Here v is the kinematic viscosity, and f = (f ,g) represents the density of body forces. 


As is well-known, the appropriate space for the pressure is L 2 (Q) defined as follows 
Lq(Q) = { <p € L 2 (Q) , J 0 «p(x) dx = 0 } . 

We denote by C” (Q) (resp. ^(Q) ) the space of all functions that are infinitely differentiable 
and are 2n-periodic in the second direction as well as their derivatives (resp. that are infinitely 
differentiable with compact support in the first direction and 2n-periodic in the second one as 
well as their derivatives). In order to take the boundary conditions (2.2X2.3) into account we 
define the spaces Hq # (Q) as the closure of # (0) into hVo), and Hq(A) as the closure of 
C~(A) into H^A). Let us define also, for any positive real number r, the space H^(0) as the 
closure of (0) into H r (Q) . 

In this framework, it has been proved in [4] that the problem (2. 1 )(2.2)(2.3) is well 
posed for any force in the dual space [H^(Q)] 2 of [Hq „(Q)] 2 , the norm of which is denoted by 
|| . ||_ 1 . More precisely we have 

Theorem 2.1 : For any f=(f,g ) in [H" 1 (Q)] 2 , problem (2. 1 )(2.2)(2.3) has a unique 
solution (u,p) in [Hq ^(Q)] 2 xL 2 (Q) , andonehas 
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(2.4) II u ||, + ||p ||<c ||f II., . 

For any 1 = (f,g) in [H*(Q)] 2 , o ^ 0 , the solution (u,p) of the problem (2. 1 )(2.2)(2.3) 

verifies 

(2.5) Hu|| a+2 + ||p|| a+1 <C|f||, . 


Let us write the dependent variables in terms of Fourier series in the periodic direction 
( we denote by i the complex square root of - 1 ) 


( 2 . 6 ) 


u(x,y) = Lm = _oo u m (x) exp(j. my) , 

P(x,y) = P m (x) exp(imy) , 

As is well known, this procedure decouples all Fourier modes for the Stokes problem , resulting in 


th 


the following set of equations for the n mode 


(2.7) , -v(u xx - n 2 u n ) + p x = ? n , 

(2.7) 2 -v(0 x n x - n 2 v") + i np n = g n 

(2.8) u" + i nv n = 0 , 


and 


| A p°(x) dx = 0 , 

with the following boundary condition on u n 

(2.9) u n (-1 ) = u n ( 1 ) = 0 . 


The following proposition is now straightforward. 

Proposition 2.1 : For any f = (f,g ) in [H (Cl)] 2 , the pair (u.p ) in [Hq # (Q)] 2 xL 2 (Q) is 
the solution of problem (2. 1 )(2.2)(2.3) if and only if its Fourier modes are solutions in 
[H^(A)] 2 xL 2 (A) for n & 0 (resp. [HQ(A)] 2 xLQ(A)for n = 0 ) of problems (2. 7) ( 2. 8) 

(2.9). 


Our numerical technique will be based on variational forms equivalent to (2.1)(2. 2)(2.3) 
and to (2.7)(2.8)(2.9), We first introduce the notations 
X = [Hj # (Q )] 2 , M = L 2 (Q) , 
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and provide X with the standard semi-norm | . |, of [ H 1 ( 0 ) ] 2 equivalent to || . ||, over X. 

The variational formulations are given by: Find (u,p) in XxM such that 
(2.10) v(Vu,Vw) - (p.divw) = (f,w) , Vw <E X , 

(2. 1 1 ) (q.divu) = 0 , V q € M , 

and ( dropping the subscript n ):Find (u,p), u = (u,v)/n [Hq(A)] 2 xL 2 (Q) such that 

(2.12) v[(u x ,w x ) + n 2 (u,w)] - (p,w x +inz) = (f,w) , V w = (w,z) € [Hq(A)] 2 , 

(2.13) (q,u x +i nv) = 0 , Vq€L 2 (A) . 

Remark 2.1 : Note that, for the case n = 0, problem (2.12X2.13) is well posed only in 
[Hq(A)] 2 xLq(A) ; in what follows we shall not consider that case for simplicity of formulation. 

2.2 The discrete formulation. 

Let K be a fixed number independent of the forthcoming parameters of discretization. We 
divide A into K subintervals A, ,...,A K , and set Q k = A k x 0. The spaces of approximation will 
consist of functions that are piecewise-polynomial over A k and trigonometric in the second 
direction. These discrete functions will be determined in order to verify problem (2. 1 0) (2. 1 1 ) in 
a discrete sense. More precisely, we shall replace the integral appearing in the L (A.) -scalar 
product by quadrature formulas associated with the Gauss and Gauss-Lobatto points. 

Let us introduce now the parameter of discretization h = (N,M) , a pair of IN 2 with N ^ 2, We 
denote by (P NK the set of all functions that are polynomial of degree less than or equal to N on each 
subinterval A k , k = 1 ,...,K ( in the case K □ 1 we simply write 1P N ). 

Next, we denote by $ M , the set of all trigonometric polynomials of degree less than or equal to M , 
i.e. 

S M = { <p(y) = LL_ m exp(imy) } ) . 

Let us define now the quadrature formulas on A and 0. We denote by (Cj , w.) for 
i = 1 ,... , N- 1 , the nodes and weights of the Gauss formula and by (£. , Qj) for i = 0,... , N , the nodes 



-6 


and weights of the Gauss-Lobatto formula. The following relations can be found in [8] for instance 
(L n denotes the Legendre polynomial of degree n) 

(2.14) V i € {1 .... ,N-1 } , = 0 and L^) = 0 ,^ = -1,^=1, 

(2. 1 5) c V <p € P 2N _ 3 , Lj = , <p(Ci)Wj = J A (p(x)dx , 

(2. 1 5) cl V <p 4 P 2N _, , Li :0 = J A <p(x)dx . 

Over each subinterval A k , k = 1 K , we then define a quadrature formula from the previous 

one's by a suitable affine transformation. Setting A k = ] a k , a k+1 [ , we define over A k , the sets 

^i,k ■ w i,k)i=i,N-i and <*tk - Pi.k )| s ojt as follows 
C f ,k = a k + ( a k + i- a k)^i + D /2 - 
(2.16) co. k = 2co./(a k+1 -a k ) , 

*ijk = a k + ( a k + i-Qk)^i + 0/2 , 

. Pi,k = 2 ei/(a k+1 -a k ) . 

The previous points lead to the two discrete scalar products on C°(A) 

(2. 1 7) g V (<p,»p) € [C°(A)] , ( ip , ) NG = L k= i Lj = i *P (^| (k )'l , (C j k ) w. <k , 

(2.17 ) CL V (ip.ip) € [C°(A)] , (*P>'p)N ( G|. = ^k=1^i=0 

Over ©, we consider the points 0j = -n + 2jn/(2M + 1 ) , j = 0 2M. They verify 

2 M 

(2.18) V ip € $ 2M , (2 tt/( 2M + 1 )) L jr0 <p(0j) = J© <p(y) dy . 

This gives us the following discrete scalar product over C°(0): 

(2.19) V (ip.ip) 6 [C°(0)] 2 , (ip ,«|i) M#e = (2n/(2M + 1)) I™,, (p(0j)ip(9j) . 

Finally we define the discrete scalar products over C°(Q) : 

( 2 . 20 ) c V (tp.ip) € [C°(Q)] 2 , (<p,<p) h<6 = (2n/(2M + 1 ))L k= i LuJ L^o <P(Ci, k , 0 j)vp(C lk , 0 j)co ik 

(2.20) eL V (<p ,f) € [C°(Q)] 2 , (ip ,'p) hGL = (2n/(2M+ 1 ))L k= i L" = q L™ o ^i.k > 0 p^i,k - 0 p5i,k 

Let us set 

x„ = xn ("' n.k ® V 2 ■ = m n <p„_ 2 , K ® s„) . 

The discrete problem is : Find (u h ,p h ) in X h xM u $uch that 

(2.21) v(Vu hl Vw) hGL -(p h ,divw) hG = (f,w) h6L , V w € X h , 

(2.22) (q,divu h ) hG = 0 , V q € M h . 
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Using the exactness of the quadrature formula (2.18) over $ M x $ M we can, as in the 
continuous case, decouple this problem into a set of (2M + 1 ) equations for the Fourier modes of 
the solution: Find (u N ,p N ) in [P N K fl H G (A)] 2 x IP N _ 2 K such that 


(2.23) 


V ^ U Nx . W xVoL +n K. W V6J-(PN. W x + i n2 V6 = ( f > W )N.6L - 


Vw=(w,z)€[[P NK nH ( !(A)] 2 


(2.24) (q,u Nx + inv N ) N _ 6 = 0 , V q € 0> N _ 2JC . 


Remark 2.2 : We state here an equivalent pointwise interpretation of these problems. Let us 
first work with the case K = 1 . We consider first a discrete problem close to (2.23) : Find u N in 
IP n (A) ft Hq(A) such that 

(2.25) (u Nx ,w x ) N6L + n 2 (u N ,w) N6L = (f,w) NGL , V w € 1P N D H«J(A) . 

We note that the products u Nx w x and u Nxx w belong to IP 2N _ 2 whence, from (2. 1 5) GL and (2. 1 7) GL , 
we can write 

(2.26) (u Nx ,w x ) NGL = (u Nx ,w K ) = -(u Nxx ,w ) □ -(u Nxx ,w) NGL . 

Let us now introduce the Lagrange interpolant Qj of the point £. , inO,...,N, i.e. the polynomial of 
IP N that verifies 

(2.27) Vr = 0,...,N , Q i (^ i .)=8. i . , 

where 8... denotes the Kronecker symbol. Taking w = Q. , i □ 1 .... ,N— 1 , in (2.25) yields, thanks to 
(2.26), 

(2.28) V i □ 1 .... ,N- 1 , -u Nxx (^) + n 2 u N (y = f(y . 

Problem (2.25) appears as a collocation approximation of the solution of the Poisson equation. 

Unfortunately, the same is not true for the Stokes problem (2.23X2.24). Indeed the 
discrete scalar product ( . , . ) NG involves the points Cj . Hence, we should introduce, as test 
function z, the Lagrange interpolant corresponding to that set of points. Such an interpolant in the 
expression [(u Nx .w x ) NGL + n 2 (u N ,w) NGL ] would not decouple the Gauss-Lobatto points. The only 
equations we can obtain, in the case K = 1 , are the following ( take w = ( Q. , Qj ) , 1 ^ i ,j < N- 1 , 
and note that p N Q ix belongs to IP 2N _ 3 so that (p N , Q ix ) N G = (p N . Q ix ) = -(p Nx , Qj) = -(p Nx , Q|) n ,6l ) 
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V i = 1 

N-l , 

, -v[u Nxx -n 2 u N ](g + p Nx (g = f(^ , 

, -v[v Nxx - n 2 v N ]Ui) + i L^d Q^iOp^LiOw., = f(^) 

V i = 1 

N-l 

V i = 1 

N-l , 

, [u Nx + i nv N ](Lj) = 0 . 


Remark 2.3 : Let us now consider the case K > 1 . Here, the interpretation of (2.25) involves the 

various virtual boundaries a k+1 = A k fl A k+1 , k = 1 K-1. Indeed taking the Lagrange 

interpolants of the points K iik different of - 1 and 1 gives 

(2.29) V k = 1 K , VI =1 N- 1 , -u Nxx U i>k ) + n 2 u N (L lk ) = f(S i<k ) , 

(2.30) V k = 1 K-1 , [ -U N xx ( a k) + n U N^ a k^ “ ^ a k^^N,k-1 + £*0,k^ = U Nx^ a k^ “ U Nx^ a k^ ’ 

[here g(x + ) (resp. g(x')) stands for lim g(t) ( resp. lim g(t) ).] 

t ^ X t ^ X 

t > X t < X 

Let us note that (2.29) is a collocation method for solving the Poisson equation, while 

(2.30) is, in a weak sense, the translation of u Nx (a k ) = u Nx (a k ) since g N k and g 0 k+1 are <3(N' 2 ). 

This condition on the derivatives means the continuity of u Nx , which is the usual condition added to 

a multidomain technique. 

For the problem (2.2 1 )(2.22) , we could derive a collocation-like interpretation of 
problem (2, 1 2)(2. 1 3) but this one is here not meaningfull on the boundary. 
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3. Error Analusis. 

The main result of this section will consist in an asymptotic expression for the error bound 
between the exact solution of (2.1 0)(2. 1 1 ) and the approximate solution of (2.2 1 ) (2.22). The 
result is optimal under the ( mild ) assumption that there exists a constant C* independent of h 
such that 

(3.1) M C* N . 

This is not a limitation for the practical cases of numerical interest, but, for the theoretical point 
of view, we give in a final remark the behaviour of the error bound we can prove in the general 
case. We don't consider the dependence of the error bound with respect to a possible growth of K 
involving a decay of the measure of the A k 's; this is now under consideration . The analysis is 
more technical since such a scheme would require a dependence of N and M with respect to the 
measure of the A k '$ and the various ratios of K, N and M would be involved in the estimate. 
Nethertheless if K is fixed or bounded, the case of dependence of N and M with respect to k can be 
handled by the same proofs as those explained in this section by using the general results of the 
appendix B. 

The analysis we are going to perform will use extensively some of the main properties of 
the Legendre basis of polynomials, Let us recall them before starting the proofs. We denote by L n 
the Legendre polynomial of degree n and recall that L n has the same parity as n and that 

(3.2) L n (-1) = (-1) n , L n (1) □ 1 . 

next we give the formulas that can be found in [8; Chapt. 2 ,§7 ] 

(3.3) (L n , L n .) = (2/(2n+ 1 )) 8 nn . , 

(3.4) (( 1 -x 2 )L n ) + n(n+ 1 )L n = 0 , 

(3.5) (n+1) L n+1 (x) = (2n+1) x L n (x) - nL n _,(x) , 

(3.6) L n (0 dC = [L n+1 (x) - L n _ 1 (x)]/(2n+l ) . 

As pointed out in [5] the existence and uniqueness of the solution of problem (2.2 1 X2.22) 
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relay on a compatibility condition between the discrete spaces M N and X N . We first confirm this 
out in subsection 3. 1 . 


3.1 Some properties of the discrete diver gence operator. 

The main result of this subsection consists in the following 
Lemma 3.1 : For any q in M h there exists a function w in X h such that 

(3.7) V k = 1 .... ,K , V i □ 1 .... ,N- 1 , div w(t ik , .) = q(t ik , .) , 

(3.8) Hwll, < C || q || . 


Proof of lemma 3,1: case K = 1. For any q in M h , we define a polynomial $ of 
[IP N fl Hq(A)] <s> S m by the following conditions 
(3.9) V i = 1 N-1 , A<I>(t; i , •) = q(Cj . •) . 

We prove in the appendix A. 1 that such a polynomial exists, is uniquely determined by (3.9) and 
satisfies the following bound 


(3.10) ||*|| 2 <C||q|| ■ 

Let us set w = (w ,z) = grad <J> , or again 


(3.1 1) w = <J> X , z = <i> y . 

As an easy consequence of (3.9)(3. 1 0) we obtain 


(3.12) 


Vial .... ,N- 1 , div w (Cj , .) □ q(tj , .) , 
II w II, < C || q || . 


Let us search now a function in (1P N <s> S M ) © [Ry such that 
(3.13) ^ g (±1,.)B-4» x (±1 f .) , M^ x (± 1 ,.) = 0 . 

To this purpose, let us write the functions $ and in terms of Fourier series in the second 
direction, we derive 

^(x.y) = Lm-_ M $ m (x) ex P(l € [P N n Hq(A) , 

\Kx,y) = Lm = _ M * m (x) exp(imy) + Xy , € IP N fl H G (A) . 

As a consequence of the fact that q belongs to Lq(Q), we obtain from (2. 1 5) G and (3.9) , 



0 = Ja dx J© dy A<J>(x,y) = 2 nf A dx ($ x ° x )(x) =$J( 1 ) -$J(- 1 ) . 


so that 

(3.14) $J(1)=$J(-1) 


Note also that due to the classical trace results ( see [ 1 4] ) we have 
(3.15) I^1)UC||«'X 2 .o<CJ<M 3/2 . ( , ■ 

Let us introduce, for any m in IN, two elements r* of CP N verifying 


(3.16) 


r;(1)= 1 , r*(-1) = 0 , r*‘C±1) = 0 


m 


r m( 1 ) = 0 .r;(-1) H 1 ,r;-(±1)=0 


and such that there exists a constant C independent of m and N 
(3.17) ( m~ 3 || r~ " || 2 + m -1 1| r" 1 || 2 + m || rl || 2 ) + (m" 3 ||r 


m 


m" II 2 + m' 1 ! r^ ' || 2 + m || r^|! 2 ) < C. 


In order to find such elements, we define, for any y in IN, the polynomials s y and s y by 
s y (x) = x 2y -( 1 /2) x 4y , 


s M (x) = [(4y+ 1 ) x 2y+1 -(2y+ 1 ) x 4y+1 ]/(4y) . 

Then it is an easy matter to check that r* = s E(m/4C ., ± s E(m/4C , ) are elements of [P N and solutions of 
the problem (3. 1 6)(3. 1 7)(remind that the constant C* was introduced in (3.1) to impose a 
relation between M and N). 


Let us define now M* by its components as follows 


(3.18) 


V m * 0 , im =$™(-1) r‘ + $™(1) r* , 
= 0 , 

x = -4,°d). 


It is an easy matter to check now that, as a consequence of (3. 1 4)(3. 1 6), (3. 1 3) is verified. Let 
us now estimate the H 2 - norm of 

We begin by the H 2 - seminorm of and compute it as follows 


IH=ll*Jf + K«lf + 11 %ir 


=L 


M 

m= -M, m*0 


( 1 /m 2 )( || r 


m 


II 2 + m 2 1| r 


II r m' ll 2+ m 4 || r; || 2 )|^ m (-1)| 5 


+ (1/m 2 )( II C || 2 + m 2 1| r;i 2 + m 4 || r; || 2 ) |$ x m ( 1 )| 2 
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=L 


M 


m= -M, m*0 


( m -3 1| r 


•• Il2 


.M 




m= -M, m^O 


+ m" 1 II I'd,' ll 2+ m II r; II 2 ) ml <*>"(- 1 )| 2 
( m -3 1| C || 2 + m - ’ll CII 2 + m||rJ| 2 )rn|$ > m (l)| 2 


From (3.17) , we derive as a consequence of [ 1 4] 

I * l 2 % < c (Zl , H m 1 4>:(-l )l 2 * Li m I $:<-) )l 2 ) 

<C(||$,(-I,. 

)H 1 /2 ,0 + 

<C||4>,II 2 , 0 , 


so that 


(3.19) IH,o<C||<J>|I 2 ,o ■ 

Using tF x (± 1 ,0) □ 0 and the Poincare-Friedrichs inequality, we deduce that 

(3.20) ||*J 1i0 <C||*|| 2 , n . 

Besides, we note that - Xy is a periodic function with zero average since is equal to 0. Hence, 
it is standard to note that 

II * - Xy || 0>n < C || - X || 0>o < C || 4> yy || 0>0 <01* | 2n , 

and from (3.1 5)(3. 19X3.20) we finally derive that 

(3.21) m 2 <C||*|| 2 • 

For the moment, we can notice that the function w defined in (3. 1 1 ) is not in X h , but only in 
[IP N ® $ M ] x [(IP N fi Hq(A)) ® S M ] . So, let us set w = w + curl 4>, we note that 
div w d div w + div (curl 40 = A<l> , 
and (3.8) is an easy consequence of (3. 1 0) and (3.2 1 ). 


Proof of lemma 3. 1 ; case K > 1 . Let us define the functions q k in L^(Q k ) as follows 

V (x,y) 6 Q k , q k (x,y) = q(x,y) - a k , a k = ( flk q(x,y) dxdy / meas(O k ) . 

Since we do not want any information for the values of q or div w on the virtual boundaries a k x 0 , 
we simply construct K functions w k on each Q k , k = 1 ,...,K , as in the previous case such that 

V k = 1 .... ,K , V i = 1 ,... ,N- 1 , divw k (t i>k ,.)=q k (t i>k ..) , 

II w k I1 1 ,o k ** c II 9k Ho,n k • 

Let us define the function w over Q as follows 

V x € Q k , w (x) = w k (x) . 
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Due to the fact that w k vanishes on the boundary of Q k , it is an easy matter to check that the 
function w satisfies 


Vk = 1 .... ,K , V i = 1 .... ,N— 1 , divwO;^,.) = q(t; ik , .) - a k . 

We only need now to add to w a piecewise linear function to find the good solution. More precisely 
we define w = (w,z) by 


V (x,y) € Q k 


k- 1 


w(x,y) = w(x,y) + a k (x - a k ) + £, s1 ex* measCAp 


[ z(x,y) = z(x,y) . 

It is an easy matter to check that w is still in Hq # (Q) since by hypothesis 
£ 

2n \ oc t meas(Ap = q(x,y) dxdy = 0 . 

The inequality (3.7) is also straightforward. 


Remark 3. 1 : We have used the hypothesis (3.1) in (3.34) and in the proof of the existence of 
functions satisfying (3. 1 6)(3. 17). In the general case, when (3. 1 ) need not hold, if we work with 
more sophisticated combinations of the L n we can verify that (3.17) and (3.20) still holds with 
CM/N 2 in place of C in (3.20) and following the same lines as in the proof of the appendix, we can 
prove (3. 1 0) with C(M/N)( 1 +CM/N 2 ) in place of C . Furthermore, it follows from [4] that 
(3.8) follows in the general case, with CM in place of C . 


3.2 Error estimate, 

Let us first put the discrete problem in an abstract formulation, in order to apply the 
standard results of Brezzi [6] concerning the approximation of saddle-point problems like 
problem (2. 1 0)(2. 1 1 )( see [ 1 0] or [5; sect. 1 ] for more details and [ 1 ; sect. 1 ] for a well-suited 
generalization for the numerical analysis of the Chebyshev spectral method ). 

Let us first define the bilinear forms 

(3.22) V (u,w) € [^(Q) 2 ] 2 , a h (u,w) = v(Vu,Vw) h6L , 

(3.23) V (q,w) € C°(Q)x ^(Q ) 2 , b h (q,w) = -(q,d1vw) h 0 . 

With these notations, problem (2.2 0(2.22) can be rewritten as follows : Find (u h ,p h ) in 
X h xM h $uch that 
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(3.24) a h (u h ,w) + b h (p h ,w) = (f,w) hCL , Vw<EX h , 

(3.25) b h (q,u h ) = 0 , V q € M h . 

The analysis of this problem will require four properties of a h and b h , that will be verified 
in the following lemmas. 

Lemma 3.2 : There exist two constants oc and y independent of h such that 

(3.26) V(u,w)€X h 2 , a h (u,w) < k|| u I1 1 II w || 1 , 

(3.27) V u<E X h , a h (u,u) > oc || u || 2 . 

Proof : The case K = 1 is standard ( see for instance [4; Prop. 1 1 1.4] ) and is based on the 
following inequalities ( see [7; Lemma 3.2] ) 

(3.28) V (<p,<p) <e 1P N x (P N , ( <p , t ) N (5L < 3 || cp || || ip || , 

( 3 . 29 ) v (p € 1P N , ( <p , <p )(g_0L ^ II *P i 

We detail the analysis of the case K > 1 . From (2.20) CL and (3.22) we have 

V (u ,w) €X h x X h , a h (u,w) = (2n/(2M + l )) L k=1 L, N =0 Vu(£ jk >9j)VwU ik - 9 j^i,k 
from (2. 1 8) we deduce 

(3.30) V (u,w) € X h x X h , a h (u,w) = j 0 dy [£ k=1 L, N =0 Vu( ^i,k -y )Vw ^i,k *y^i,k ] • 

Using (3.28) we obtain now 

V (u ,w) € X h x X h , a h (u,w) < 3 j 0 dy [£ k=1 1 Vu(.,y) || Ak H Vw(. ,y) || Ak ] , 

< (3/2)J e dy [L k K =1 || Vu(.,y) ||^] 1 /2 [ II Vw(.,y) igj 1 ' 2 , 

< (3/2) <J e || Vu(.,y)|| 2 dy) 1/2 (} 0 || Vw(.,y)|| 2 dy) 1/2 , 

<(3/2) Hull, || w ||, , 

which proves (3.26). From (3.29) and (3.30) we obtain (3.27) by using the same arguments. 

Let us now analyze the properties of the discrete bilinear form b h . 


Lemma 3.3 : There exist two constants g, 8, independent of h such that 

(3.31) V (q,w) € M h x X h , b h (q,w) < 8 || q || || w ||, , 
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(3.32) Vq€M h , sup b h (q,w) > & || q || || w||, . 
w€X h 


Proof : From (2. 1 8)(2.20) G and (3.23) we deduce 

K N - 1 

V (q,w) € M h x X h , b h (q,w) = J 0 dg [£ ka1 L i= 'i q(i; lk .y) divw(t jk ,y)o) ik ] 
Using the Cauchy-Schwarz inequality gives for any (q,w) in M h x X h 


,N -1 


b h (q,w) < } 0 dy L k=1 ([£," ~i q 2 (tj <k .y)u iik ] 1/2 [Lj7i (div w ) 2 (Ci, k ,y)co ik ] 1/2 ) , 


,N -1 


,1/2n 


,N -1 


N -1 


< [Je dy Zl k= 1 Z^i q (^ k ,y)o i>k r 2 [Je dyL^L^", (divw) 2 (t jk ,y)u. k ] 1/2 ), 


Let us notice that q 2 (.,y) is in IP 2N _ 4 ; (2. 1 5) 6 yields 

(3.33) b h (q,w) < ||q|| [J 0 dy L k=1 (divw) 2 (t; ik ,y)co ik ] 1/2 , 

Writing now divw in the Legendre basis over A k , k = 1 ,...,K 

(3.34) V x € A k , V y € 0, k = 1 K, divw(x,y) = =0 i nk (y) L n (x) - 


gives 


,N -1 


N-1 


£ k = 1 2 -i = 1 (divw) ( C j k . y ) to j k - ^ kl i Ljj 1 (ZL n -o T n,k(y^ *- n (t jk ) ^ W i,k 


.N-1 


and finally 

(3.35) L k=1 Lu 1 1 (divw) 2 (C; ijk ,y) w ik = L k=1 Ln =0 Llo( T n> T vA^ y ^^i = i ( L n L v)(^,k) w i. k ] 
As a simple consequence of (2. 1 5) 0 and (3.3) we deduce 

(3.36) V (n,v), n + v < 2N-3 , (L n L v )(q ik )co. k = (L n , L v ) = 8 nv (2/(2n+ 1)) 

Thanks to (2. 1 4) and (3.5), we derive that 


V i = , L N (C i/k ) = -((N-O/N) L N _ 2 (C lk ) , 

V i = 1 .... ,N- 1 , L N _ 1 (C ik ) = 0 , 
whence, using again (2. 1 5) 0 gives 

<l n l n ) <c u ) =«n-o/n ) 2 i";,' (l„. 2 ) 2 (Ci*) o lk 

= ((N-1 )/N) 2 J a (L n _ 2 ) 2 (x) dx 
< 4/(2N+ 1 ) = 2 J a (L n ) 2 (x) dx , 

£-i = 1 (l-N L N-1^(((i,k^ W i,k = 0 ’ 
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. N - 1 


,N -1 


= 1 ^ L N L N-2^^i,k^ 10 i.k - ~((N-1)/N) 2Li = i 0- N-2 ^ U i ( k 

= - ((N-O/N) J a (L n _ 2 ) 2 (x) dx 
< 2 J a (L n ) 2 (x) dx , 

= 1 ^*- N _ 1 L N _ 1 )(^ !<k )(0 j(k = 0 . 

The four previous equalities together with (3.35) and (3.36) yield 


.N -1 


.K _ N 


Li = , (divw) (C iik ,g)W lk < 3 L k = 1 L n -_ 0 L v= o [JA k < L n L V >< X > dX ]> 

3 f A (divw) 2 (x,y) dx , 

which, from (3.33) gives (3.3 1 ) with 8 = V3. 

Let us now prove that the compatibility condition (or inf-sup condition) between the 
discrete spaces M h and X h is satisfied. From Lemma 3, 1 we know that there exists a function w in 
X h such that 

(3.37) V k = 1 K , V i = 1 .... ,N— 1 , div w(t i<k , .) = q(C ik , .) , 

and 


(3.38) || w||, < C H q II • 

With (2. 1 5) 6 and (3.37) we easily verify that 
b h (q,w) = (q,q) N G = || q|| 2 , 
using now (3.38) yields 

b h (q,w) > C || q || | w H, . 


The previous lemmas prove that the approximation of the Stokes problem by the scheme 
(2.2 1 X2.22) is well-posed. More precisely we obtain 

Theorem 3. 1 : For any f in (C°(Q)) 2 , there exists a unique solution (u h , p h ) to problem 
(2.20(2.22). Moreover, if we assume that hypothesis (3.1) holds and that f belongs to 
H*(Q) 2 , a > 1 , the following error estimate for the velocity and the pressure holds for any 
y > 1/2 

(3.39) || u - u h ||, + || p - p h || < C (M- ff + ( 1 + M- y N V /2 '0 || f II, . 
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Proof : The existence and uniqueness of u h and p h follows directly from Lemma 3.3 and 3.4 and 
[10, Theorem 1.1]. Moreover another consequence of that theorem is the following 

(3.40) ||u-u h || 1+ ||p-p h ||<C[ inf ||u-w|| 1+ inf |[ p - q|| + sup h 

w€X h q€M h w€X h 1 

and the result is an easy consequence of Theorem 2. 1 , Theorem B.5, Corollary B. 1 and Theorem 
B.8 of the Appendix, 


Remark 3.2 : Note that Theorem 3. 1 still holds when f only satisfies 
f| Qk e H“(n k ) , 
see appendix B. 
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4. NUMERICAL RESULTS. 

4.1 Implementation . 

In this section we describe the details of the discrete equations (2.23) and (2.24). We start 
by defining the bases for the space [IP NK Pi H G (A)] 2 x 1P N _ 2 K in which we search for our solution 

(u N ,p N ) . As described in Section 2.2 A is divided into K spectral elements A, A K .lneach 

element A k the velocity u N of [1P N K f| H G (A)] 2 is expanded in terms of N lh order Lagrangian 
interpolants Q. (see (2.27)) through the Legendre Gauss-Lobatto points . We then define a 
mapping § from x € A k onto r€l=]-1,1[asr = -1 + 2(x - a k )/(a k+1 - a k ). Then we state 

(4.1) V x € A k , uj(r,.) = u N (x,.) 
and 

(4.2) u^(r,.) □ £, N _ 0 uj^Q^r). 

Here u k = u N (£ jk ,.) is the velocity at the (local) point £ ik in the interval A k ; that is, 
(4. 1 )(4,2) is a nodal basis. Similarly to the velocity the data f is also expanded in terms of N th 
order Lagrangian interpolants through the Legendre Gauss-Lobatto points , 

(4.3) fj(r) = L" z0 ffQi(r) . 

The pressure p N € IP N _ 2K is expanded in terms of (N-2) lh order Lagrangian interpolants Q. 
through the Legendre Gauss points Cj . 

(4.4) pj(r) □ Pj k Q, 

where p k = p N (t jk ) is pressure at the (local) point t; ik in the interval A k . Note that the Gauss 
points are naturally suited for the pressure, which need not be continuous across elemental 
boundaries. 

The expansions (4. 1 )-(4.3) are now inserted into (2.23) and (2.24) , and the discrete 
equations are generated by choosing test functions w € [IP NK Pi Hq(A)] 2 in (2.23) which are 
unity at a single C ik and zero at all other Legendre Gauss-Lobatto points, and test functions 
q € tP N _ 2 K in (2.24) which are unity at a single £ ik and zero at all other Legendre Gauss points. 
To evaluate the integrals in (2.23) and (2.24) we use numerical quadrature through the Legendre 
Gauss-Lobatto points £ ik and the Gauss points C ik . denoted ( . , . ) N GL and ( . , . ) N c 
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respectively. The first term in (2.23) can then be written as 

(4.5) (u Nx ,W x ) N<(3L = L k= K , (L k /2) L- =0 Ho Pm D mi D mj u f ■ 
where the derivative matrix D is defined as 

(4.6) Dp, = (dQ,/dr)(g. 

We recall that p m are the quadrature weights associated with the Legendre 6auss-Lobatto points 

I 

£.,al$oL k =(a k+1 -a k ) and£ denotes direct stiffness summation. 

The second term in (2.23) becomes 

(4.7) ( u n> w Vgl = ^k = i (*-k /2)pju|‘ , 

while the right-hand side of (2.23) can be written as 

(4.8) w Vgl = £k = i (L k /2)Q.f! c . 

The left-hand side of (2.24) becomes 

(4.9) (q,u Nx + inv N ) N0 = (o. (Lj D.j uf + in L k /2 iy vj) 

where the derivative matrix t> and interpolation operator f are defined as 

(4.10) t> pq = dQ q /dr (Cp) , 

(4.11) r„, = <yg. 

Note that in (4.6) no direct stiffness summation need be performed since the Legendre Gauss 
points t; ik are all distinct. 

In matrix form the set of discrete equations (2.23) and (2,24) can be written as 

< 4 - ,2 > a w u - g .p P p = bi 
( 4.13) D w u=0 

where A,,, is given bg (4.5)-(4.7) , B lpp bg (4.8) , D <pp bg (4.9) .and G app is the 
adjoint-matrix of . The Uzawa method used to solve (4. 1 2) -(4. 1 3) will be described in more 

®PP 

details in a future paper[ 1 9]. Basically, it consists in solving the following zero lh -order equation 
for the pressure 

D .PP (V' G .ppP = °,pp (A,pp)'' G f . 

by using a conjugate gradient algorithm , and then recovering the velocity from (4. 1 2). 
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In this section we discuss some numerical results obtained by solving the set of discrete 
equations (4. 1 2)(4. 1 3). Two test problems have been solved, both in which v = 1 and n = 1 . 
The solution (u,p) and the data f = (f.g) for the first test problem are 
f u = - ( 1 + cos n x)/n , 


(4.14) 


v = i si n n x 


sin n x 


f = - ( 1 + cos Tt x)/n , 
g = i(2+n 2 ) sin n x , 

while the solution and the data for the second test problem are 
f u = -( 1 + cos n x)/n , 


v = i sin n x , 

(4.15) p = sin n x | x— 1 /2 | v+2/3 , 

f = [- ( 1 + cos n x)/n + y+ 2/3) $gn(x- 1 /2) | x— 1 /2 | v+2/3 ] , 
g = i [(2+n 2 ) sin tt x + | x- 1 /2 | v+2/3 ] . 

Note that the solution and the data in the first test problem are infinitely smooth, while the 
regularity of second test problem is determined by the value of * , which is assumed to be an 
integer. 

In the first test problem A is divided into 2 equal subintervals A, and A 2 , i.e. k = 2 , 
while in the second test problem only one element is used, i.e. K = 1 . The numerical solutions are 
compared with the analytical solutions for different values of N , the order of the polynomial 
expansions (4. 0(4.3). To measure the error in the numerical solutions, the following error 


measures are used 


(4.16) ||u-u N || I0L = {L k K s1 (L k /2) [L t "=o [utt lk )-uf] 2 )]} ,/2 

(4.17) ||p-p N || M = (Z: k K =1 (L k /2)(L 1 N :, , a.,[p(C i , k )-p, l ] :! ))' /2 . 


Ai2m1/2 


In the first test problem we obtain exponential convergence as the order N of the polynomial 
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expensions is increased. Figure 1 shows the error in the velocity and the pressure as a function of 
the total number of degrees of freedom (Legendre Gauss Lobatto points) in the x-direction. The 
rapid convergence rate is expected due to the fact that the solution is analytic. 

In the second test problem we obtain algebraic convergence as the order N of the polynomial 
expansions is increased. Figure 2 shows the error in the pressure as a function of the total 
number of degrees of freedom (Legendre Gauss Lobatto points) in the x-direction for kb 3 and 
K = 5 . The convergence rate is given approximately as N' (v+1) . Although the error estimates 
(3.40) is somewhat pessimistic as regards the error due to the forcing term (f € H v => 
Up— p N || < N 1_v ) , as regards the approximation errors (p € H v+1 => ||p-q N || < N" 1_v ) the bound is 
quite tight. 
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APPENPIX, 

Appendix A - Approximation o f a Discrete Lapla ce Equation. 

The proof of the compatibility condition between the spaces of velocity and pressure 
involves some results concerning the approximation of the solution of the Laplace equation by a 
collocation method based on the Gauss points. 

Lemma A. 1 : For any q in M h there exists a unique <J > in [P N D Hq( A)] ® S H $uch that 

(A. 1 ) V i = 1 N-1 , ASttp.) =qtt,,.) • 

Proof : Let us consider the collocation problem : Find $ in [IP N fl H^( A ) ] ® S M such that 

(A. 2) V 1 = 1 N-1 , V j = 0 2M , A^(Ci , 0j ) □ q(C t , 0j ). 

It is an easy matter to check that $ satisfies (A.1). Multiplying both sides of (A. 2) by 
(2n/(2M + 1 )) ^(Cj , 0j )u. and summing up with respect to i and j leads to the equation 
(A. 3) V*€(P N ©S M , (A4>^) h0 = (q^) h(3 . 

Taking equal to the Lagrange interpolant of the point (C ( , 0j ) in [IP N fl Hq(A)] © S m , proves 
that the problem (A.2) is equivalent to: Find $ in [ IP N H Hq ( A ) ] <s> S M such that 
(A. 4) [IP N D H^(A)] <s> S M , (A<?^) h6n (q^) hG . 

Let us set 

(A. 5) V ('±',X)€ {[IP n PI H^(A)] © S M } 2 , c(V. X) = -(AV, X) h<(5 

In order to prove that c is continuous and elliptic, we recall that (see [2; Lemma II 1. 1 ]) 

(A. 6) V ip 4 IP N _ 2 . J A (1-x 2 ) 4r(x) 2 dx < (( 1 -x 2 ) ip ,4») N G < 2 f A (1-x 2 ) ^(x) 2 dx , 
and that (see [2; Lemma III. 2]) 

(A. 7) V * € IP N . 2 , C N -1 J a ( 1 -x 2 )%(x) 2 dx < (( 1 -x 2 ) 2 * ,4<) N 3 ^ C J A ( 1 -x 2 ) 2 *(x) 2 dx . 

For any ^ andX in [IP N fl H G (A)] © S M we derive from (2.1 5) G and (2. 1 8) that 
(A.8) c(*,X) = J 0 [-(V xx (. I y) l X(.,y)) Nf6 + (V u (. l y),X u (. > y)) N<e ]dy. 
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Let us write ^(..y) and X (,,y) in the basis ( 1 -x 2 ) L n , n = 1 

V(x,«) = L„ n :,’ t n (g) (i-x 2 ) L(x) , 

X(x,y) = LnI 1 1 X n (y) (1-x 2 ) L n (x) . 

Using (3.4) gives 

-0±' xx (.,y),x(.,y)) NG = Lnl! L"l\ v n (y) x,(y) <(< + 0 (d-x 2 )L n ,Lj ) N(5 . 

From (2. 1 5) 0 (3.3)(3.4) and (A. 6) we deduce 

(A. 9) -(^ xx (.,y),X(.,y)) N6 < (4/(2n+ 1 )) Y„(y)x„(y) (n(n+ 1 )) 2 < 2 } A ^ x X x (x,y) dx 

and that 

(A. 1 0) -(^ xx (.,y),'i'(. > y)) Ne ^ (2/(2n+1 ))£„"] <V n (y) 2 (n(n+ 1 )) 2 ^ -J A * XX (x,y)¥(x,y) dx 

Let us write now 

M/yCx ,y) = ( 1 -x 2 ) *(x,y) , 

X g (x,y) = (1-x 2 ) ^(x,y) , 
using (A.7) yields to 


(A. 1 1 ) (%(.,y),X u (.,y)) N0 <C (* y (.,y),X y (.,y)) , 

and to 

(A. 1 2) (M/yC. ,y) ,M/ y (. ,y)) N G > C N" 1 (¥ g (. ,y),^ g (. > y)) , 

Finally, due to (A. 8) (A. 9 ) (A. 1 1) and the Poincare-Friedrichs inequality, we deduce that c is 


Ml I I I Ml II I 


]y continuous over [!P N f] H q ( A ) ] ® S M . Due to (A.8)(A. 1 0)(A. 1 2) , we derive a 
(nonuniform) ellipticity of c over [IP N fl H G (A)] <s> $ M , more precisely, we obtain 
(A. 13) [P N nH’(A)]< 3 )S M , c(¥ 1 ¥)>CN“ 1 |M|? • 

From the Lax-Milgram lemma we conclude that problem (A.4) is well-posed, i.e. (A. 4) admits a 
unique solution $ in [IP N fl H G (A)] ® S M . Besides from (A. 1 3) we derive the following estimate 
(A. 1 4) 1*11, <CN- 1 ||q|| . 


In this second lemma, we are going to derive a uniform bound for the H 1 (0)-norm of $ x . 
Lemma A. 2 : The solution * in [IP N fl Hg( A)] ® S M of problem (A. 1 ) satisfies the following 
estimate 

(A. 15) ||* x || 1 <C||q||. 
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Proof : Let us choose now ^ = <P XX in (A. 3), we derive 

^xx ’^xx^h.G + ^yy’^xx^h.G - ^’^xx^h.G ’ 

and we obtain from (2. 1 5) G , (2. 1 9) and (2.20) 0 that 
(A. 16) ll<!>„l| 2 -( 9 (<I> u .1>„ s ) N , 6 <llq||l|f >< ll ■ 

Writing^yC-.g) in the basis ( 1 — x 2 ) , n = 1 , 

(A. 17) V (x,g) € 0, 4> y (x,g) = Li"] <P n (y) ( 1 -x 2 ) L n (x) , <p n € S M 
gives, as in (A. 1 0), 

(A. 18) -(‘Py/I’xxyVo ^ (2/(2n+1)) LH;1 4> n (y) 2 (n(n+1)) 2 ^ J A <P x 2 y (x,g) dx . 

Finallg, we conclude with (A. 1 6) that 
(A.19) II^H, < 2 || q || . 

Unfortunatelg, the inequalitg || $yy || < C || q || is not so easg to derive. This is done is the 
following lemma. 

Lemma A, 3 : The solution $ in [1P N fl H G (A)]<*> $ m of problem (A. 1) satisfies the following 

estimate 

(A. 20) \\* v l <C||q||. 

Proof : We first define a function such that 
(A. 21) V 1= 1.....N-1 , $(£.,.) = <Kt,,.) , 

and such that the inequalitg || $ yy || < C || q || holds. To this end we use the basis ( 1 — x 2 ) L n , 
n = 2.....N, 

(A. 22 ) V(x,g)€0, 4>(x,g) = i: n N =2 T3 n (y)0-x 2 )C(x) , T, n € $ N , 

From (3.4), (3.5) and (3.6) we obtain that 

(A. 23) f (1— x 2 ) L n '(x) = (1-x 2 ) L n '_ 2 (x) - N(2N-1) L N (x) + N(2N— 1) x L N _,(x) , 

( 1 -x 2 ) L n _, (x) = 2 x L n _,(x) + N(N- 1 ) L N-1 (x) , 

so that 

V i = 1 .... ,N— 1 , ( 1 — Cf)L^(C i ) = ( 1 -x 2 )L N _ 2 (t: i ) - N(2N-1 )L n (Cj) , 
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V 1 = 1 N-1 , (1-Cf)L N '. 1 (C i ) = 2C i L N _ 1 (i; i ) • 

Hence the function $ defined by 

(A. 24) V (x,y) € Q, $(x,y) = i] n (y) ( I -x 2 )L n '(x) + 2Tj N _ t (y) xL^/x) 

+ti n ( y)[( 1 -x 2 )L n _ 2 (x) -N(2N-1)L n (x)] , 


satisfies (A.2 1 ). 

Moreover we have 

(A. 25) V (x,y) € 0, $(x,y) = < ?(x,y) - ti n _ 1 (y) (N(N-1 )) L N _ t (x) 

- %(y) (N(2N- 1 )) x L n _ 1 (x) 


Using now (A. 3) and (A.2 1 ) gives 

^xx ’^gy^h.G + ^yy’^yy^h/G ~ ^’^yy^h.G ’ 

so that as previously from (A. 1 8) 

(A. 26) IH> >S II 2 + ( VwVc < II q II 2 ■ 

From (A. 24), we obtain 

(A. 27) $ yy (x,y) = Vu) d-x 2 ) L n (x) + [(Wy) + Vy)Hl-x 2 ) L n _ 2 (x) 

+2 73n_i ( y) x L^/x) - t) n (y) (N(2N- 1 )) L N (x) . 
Using now the formula derived from (3.3) , (3.4) and (3.6) yields for any n ^ v 
J A (1-x 2 ) 2 L n '(x)L^(x) dx = J a (2x L n (x) - (n(n+1)) L n (x))(1-x 2 ) L v (x) dx 

= f A [(( 1 -x 2 ) L’ (x)) (2 xLjx)) - (n(n+0) L n (x) (1-x 2 ) L v (x) ] dx 
= J A - (n(n+ 1 )/(2n+ 1 )) (L n+1 (x) - L^/x)) (2 x L v ’(x)) dx 

- f A (n(n+ 1 )) L n (x) ( 1 -x 2 ) l”(x) dx 
= J A (n(n+1)/(2n+0) (L^x)) (2xL^(x)) dx 

- | A (n(n+ 1 )) L n (x) (1-x 2 ) L v (x) dx 


To derive the value of X n we only calculate the leading coefficient of 2 xL n (x) and ( 1 -x 2 ) L n (x) , 
this is done by using (3.5) , we obtain 

\ n = J A (n(n+ 1 )/(2n+ 1 )) (2(2n- 1 )(n- 1 )L n . 1 (x) + ) 

+ (n(n+1)) L n (x) (n(n-1) L n (x) + ...) dx 


and by (3.3) 
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X n = 2 (n-1)n(n+1 )(n+2)/(2n+1 ) . 

The case n < v is treated in a symetric way. We can write, for any n and v in IN 
J A 0- x2 ) 2 L>)L>)dx = X n 8 nv . 

We compute from (A.23) and (A.27) on the one hand 

(A. 28) || $ yy || 2 = L n _2 II T) n II 2 X n + II TIn_2 + H 2 ^n-2 + 4 II "On_i II 2 II X L N _ 1 l| 2 

+ || i] N ' || 2 2[N(2N- 1 )] 2 ( 1 /(2N+ 1 )) , 

and on the other hand from (2. 1 5) 6 

(A. 29) ^yy’^yy^.G - ^n = 2 H H + II ^N-2 + II ^N-2 + ^ II T)|\|_i II ( x ^- N _ 1 , L^_ j ) NG 

+ 11 T) n ' II 2 [N(2N-1)] 2 (L n ,L n ) ng 

2 (N(2N- 1 )) ( T)n _2 + T)n >T)n )<=>((! — x ) Lfj _2 , )^ G . 

From (2. 1 5) G , (3.3) and (3.5) we deduce 

(A.30) (Ln.LnV^O-N- 1 ) 2 (L N _ 2 ,L N _ 2 ) N6 = d-N' 1 ) 2 (L n _ 2 ,L n _ 2 ) ^ 1/(2N+1) . 

Using now (3.6) combined with (3.5), we compute 

xL b ., = (N-l) l„_, * (2N-5) L n _, +CN-4) L N . 5 * xL N . 5 , 

so that 


< x2l n-i .L N -i) N .g = <2N-5) 2 [2/ (2N— 5) ] + ||(N-4) L N _ 5 + xL N _ 5 || 2 


and 


II xL n _, || 2 - (N-1 ) 2 [2/(2N- 1 )] + (2N-5) 2 [2/(2N-5)] + ||(N-4) L N _ 5 + xL N _ 5 || 2 . 
We conclude with 

(A. 3 1 ) || xL N- 1 || < 2 (x L N-1 , i )m Q • 

Finally, from (3.5) and (2. 1 5) 6 , we derive that 

(( 1 -x ) L n _ 2 , l_ N ) N 6 = ((N-l )/N) (( 1 -x 2 ) l _ N _2 , L n _ 2 )|y| e 

= ((N— 1 )/N) (( 1 -x 2 ) L n '_ 2 ,L n _ 2 ) . 

Next, using (3.3), we obtain from the equality ( 1 -x 2 ) L N '_ 2 = -(N-2)(N-3) L N _ 2 + .... 

(A. 32) (( 1 -x 2 ) L n '_ 2 , L n ) N6 = - 2(N-2)(N-3)(N- 1 )/N(2N-3) = O(N) 

Recalling now (A.26)(A.28)-(A.30), we prove that 

(A. 33) II + yy II < C ( || q | + N 3/2 II D „. 2 * D N ' II* N 3/2 || % II ) 

Formula (3.6) combined with (A.1 7) and (A.22) yields 
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(A.34) | T) N _ 2 (y) + Tj N (g) = tp N _ 3 ( y) /( 2N-5) , 

TD N -i<y) = <P N _ 2 (y)/(2N-3) , 
ij N (y) = <P N _,(y)/(2N-0 . 

In order to estimate ||<P yy || we note that, from (3.3)(3.6) and (A.25) we obtain 


II%II«c[||%|| + n 3/2 || -o N ., || * n 3 ' 2 || i)j|] 


, 3/2 


and from (A.33)(A.34) we deduce 

ll%NCtl|q|| + N ,/2 [|U„_, 


II <Pn- 2 II + II fN-1 HI 


Let us now use the inverse inequality 

V (r,s) € \R 2 + , 0 < r «; s , Vx«S M , || X ll s < C M s ' r || x || r , 
it follows that 


11% II < c III q II + M( N 1/2 [ II || ♦ 1| II + II 111 )] 

It is an easy consequence of (A, 1 7) that 

<f> n (y) 2 (n(n+1)) 2 2/(2n+1) ^ J A <P x 2 y (x,y) dx , 
so that, from (A. 1 9) we derive 

N 3 ' 2 [ IU N _, II ♦ l *„_ 2 IMI ♦„.,«]< Oil q II , 


hence, we conclude 


(A. 35) l|* gy ||<C(1+M/N)||q|| , 

which, thanks to (3.1) and (A. 1 9) , finishes the proof of Lemma A. 3 , 


As a conclusion, we state 

Theorem A. 1 : For any q in M h there exi$t$ a unique $ in [IP N f| H<J(A)] <*> S^such that 

Vi= 1 N-1 , A^CCj , .) = q(Cj , •) . 

Moreover, satisfies the following bound 
(A. 36) ||*|| 2 <C||q||. 
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Appendix B- Error Bound for the Projection and Interpolation 
Operator s 

The final error estimates require some technical results about the orthogonal 
projection operator for the Hq ^(Q)-scalar product onto the space 1P NK <s> S M and interpolation 
operators with values in the same space. 

B.1 The one-dimensional case. 

In this paragraph, we extend the results of [7] and [16] to state some properties of the 
approximation operators over P N K . See also [9] for some partial results in this direction. 

The possibility of using different values of the parameter N in each subdomain Q k was only 
evoqued in the previous part of the numerical analysis. The only difficulty that this should imply, 
would be to complicate the reading of the proofs, and absolutely not of mathematical nature. 
Howether , here we shall consider such an eventuality since the extension from the case where N is 
assumed to be constant, to the one where N is variable, is not straightforward. The interest of 
doing so is to be able to fit the regularity of the solution and, in particular, to increase the number 
of degree in the region where the solution is a bit less regular. This is a first step toward the 
general situation; the second one will take into account the size of A k and the possibility of taking 
the parameter K as a discretization parameter. 

Let us define for each k, 1 < k < K, an integer N(k), that will take now the place of the 
previous notion of N, the degree of the polynomial in the nonperiodic direction; the corresponding 
space of polynomials over A will be noted 0P N K (A) and will consist in 
= { 'P € 1- (A)> f |A k € ^N(k)(^k) } 

(note that from now on , we shall precise the interval where the variable are defined for the 
various spaces of polynomials). 

The regularity of the solution being possibly different on each A k , we introduce some spaces with 
broken norms. More precisely, for any K-tuple of positive real numbers rn (r 1f r 2 ,...,r K ), we 
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define 

H r (A) = { <p € L 2 (A), <p, Ak € H k (A k ) } , 
this space being equipped with the norm 

V ip € H r (A) , IUII r . A = [L k K ,, IUII r \. Ak l ,/2 . 

First of all, let us consider the L 2 (A)-projection operator TT n onto (P NK (A) . We have 
Theorem B. 1 : Let r be a K-tuple of nonnegative real numbers. We have for any function ^ in 
H r ( A) 

(B.1) Wv - TT N ^|| 0 A <C[Z : k K =1 N(k) k 11^ ll? k , Ak ] 1/2 • 

Proof: First, we note that from the definition of Tt N , we have 
(B.2) V^€L 2 (A) , V<P<eP nk (A) , (^-TT n ^,<J>) a = 0 , 

so that for any k in DM , 1 < k < K, we have 

V 'i' € L 2 (A) , V 9 t P N(kl (A k ) , J At C <i'-TT N <i«)| Ak <t> dx = 0 , 

hence we note that ( TT N M/ )| Ak is the projection of ^| Ak onto respect to the 

L 2 ( A k )-scalar product. As a consequence of classical results (see [7; Theorem 2.3]), we derive 
that for any r k > 0 , we have 

V * € H\A k ) , || 't'-TT^ \\ 0Ak < C N(k) k || ¥ || rk>Ak . 

Summing up the square of these inequalities, we deduce that 

v*€H r (A) , ll^-TT N ^|| 2 A = L k=1 ||^-Tr N ^|| 2 Ak 

^Cl k K =1 N(k)' 2rk ||^|| 2 k>Ak , 

and (B.1) is proved. 


Next, we state the following inverse inequality 
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Lemma B . 1 : Let rand s be two K - tuples of real numbers, such that for any k, 1 < k < K, 
0 < r k < s k . We have for any any function V in tfy K <A) 

4( s p ) 

(B.3) II 4- |f, A s C N(k) k k 


2(cr— q) 


W 


•eJ a,b[ 1 


Proof : Here again this result is a simple consequence of the following classical inverse inequality 
over (P N (]a,b[) for any a and b in IR : 

(B.4) V w 4 (P N (]a,b[) , V(q,c) € IRj, q s? o , || w || ff ]ab[ < C N 
Indeed, we have 

Vv«V<A) . Ilv ||* A = Lk., II V Ilf t Ak . 

from (B.4) applied on each subinterval A k , we deduce that 


2 


v v € ^n,k ( a ) > IMI 9 2 , A <cL k=1 N 
and (B.3) is proved. 


" r k’ A k ’ 


Now, we are interested in some projection operator from Hq(A) onto CP N K (A) fl H<j(A). 
Theorem B.2 : There exists an operator TT^ from Hq(A) onto ff> NK (A) fl H<J(A) verifying 
for any function 4* in R 9 (A) fl Hq(A) , with s being one K-tuples of real numbers ^ 1 
(B.5) V r = ( r k )i ^ k ^ K > ® ^ ^ ^ > II ^ - Hr,A ^ ^ [£-^=1 ^00 


2 ( V s k ) 




1/2 


Proof : Let us recall that, for any a and b in IR, there exists a projection operatorn N from 
H 1 ( ]a,b[ ) onto IP N (]a,b[) satisfying (see [3; Corollary IV. 2]) for any 0 < r < 1 < s 
(B.6) V w € H s (]a,b[) , || w - n N w || r ]ab[ < C N r ‘ s || w || s ]ab[ , 
and 

(B.7) n N w(±1)ow(±1). 

Let us define the projection operators n N(k)k , for any k in DM, 1 k < K, as being the projection 
operators from H’(A k ) onto IP N ( k )(A k ). We deduce from (B.7) that the element TT^ ^ defined on 
each A k by 
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VX€ A k , TtJ'Kx) = n N(kU (^ |Ak )(x) 
is an element of [P N(k) (A k ) fl Hq(A) that satisfies, due to (B .6) 

\W- Tti V ll r , A < c [L k K =1 N(k) 2<VSk) 11^ || s 2 k;Ak ] 1/2 . 

Then, (B.5) is proved . 

Let us analyze now some properties of the operator of interpolation i N K in P N K ( A) over the 
Gauss-Lobatto points. Since the degree N(k) of the polynomials of approximation can vary with k, 
we must redefine the points of interpolation. It consists over each A k of the (£jk)j : o,N(k) defined in a 

similar way as in (2.16). Using the same techniques of decomposition of the interval A in 



U k _i A k , we deduce from the classical results on the operator of interpolation in P N (A) over the 
Gauss-Lobatto points (£;.), i=0,..,,N ( see [7; Thm. 3.2]) that 

Theorem B.3 : Let r be a K -tuple of real numbers, such that r k > 1/2 .We have for any 
function ^ in Fl r ( A ) 

(B.8) ll'4’-i N *'l'll 0 , A SC[L k K , 1 Wk)’' 2 ' 1 imi r 2 l , Ak ]' /2 . 

B.2 The two-dimensional case. 

In this paragraph, we shall combine the results of section B.1 with the classical results 
concerning the approximation theory related to the Fourier case. These results and the techniques 
we use are very close to those of [4; Appendix] and [17], 

As in the previous section, we shall consider that the regularity of the function we want to 
approximate is different on the various Q k . To this hand, we associate with each Q k a couple 
(N(k), M(k)) of integers and consider the space of approximation 

<^h = { <P € L 2 (Q) , V k, 1 < k < K , <P |o k € ^NdO^k^ ® ^M(k) ) • 

Then, for any K-tuples r and s of positive real numbers, we consider the spaces 

H r,s (Q) = { <p € L 2 (Q), V k, 0 =$ k < K, <p , 0k € H rk (A k ;L 2 (©)) fl L 2 (A k ;H* k (0)) } . 

As in section 2 , we define also the spaces R r "*(Q) as being the closure of C“(Q) in R r ' s (Q). We 



-32- 


shall use in the proofs some norms over Q k , the space H r (A k ;L 2 (0)) is provided with the norm 
|| . Il r0 0k ; the space L 2 (A k ;H*(0)) is provided with the norm || || 0 s0k . 

First of all, let us consider the L 2 (Q)-projection operator TT h onto <s& h in L 2 (Q). We have 
Theorem B.4 : Let r and s be two K-tuples of real numbers. For any function tF in H r, *(Q) 
(B.9) || * - Tty || 0 . A S C [ZL (NOc)' 2 " 1 II V | r 2 0A + M(k)‘ 2!k || ¥ || 0 2 SkA ] m 


Proof : Let us denote by TT^ the L 2 (0)-projection operator onto S M (0). We recall the following 
inequality, valid for any a ^ 0 ( see [24]) 

( B . 1 0) Vv€H *(©), ||v-TT^ v|| oe < C M" s ||v|| se . 

Then, as in the proof of theorem B.l , we derive from the definition of the L 2 (Q) -projection 
operator that TT h coincides over Q k to the standard projection from L 2 (Q k ) onto 
P N / k )(A k ) ® S M(k) . It was proved in [4; theorem A. 1 ] that 


-2r. 


-2s. 


II * - Pk * lo.o k « C [N(k) I1 1- H 2 lA „ k * M(k) k II * ll 0 2 , !k , Ak ] 
we deduce that 


1/2 


ll'4'-"n h 'i'llo,o«L t= ,ll'i>-|s k >Hi 2 


-2r, 


-2s. 


<C[L l=l (N(k) k II 'i' llr k , 0 , 0 k + kick) II V IIS >Sk , Ak )l 


1/2 


and (B.9) is proved. 


In the case where N and M are constant, we have simply for the operator TT h of projection 
over P N (A) ® $ M . 

Corollary B.l : Let r a nonnegative real number. We have for any function ^ in H^Q) 

(B.l 1 ) ||^-rr h S^ I| 00 < C (N" r + M- r ) M r ,a . 

Proof : It suffices to notice that H^(Q) coincides with H r (A; L 2 (0)) f] L 2 (A ; H^(0)) (see [ 14; 
Chap. 4, Proposition 2.3]). 


Next, we state the following inverse inequality. We first define, for any K-tuples r and s of 
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positive real numbers, the space 

2e r ' s (0) = { <p € L 2 (Q), H rk (A k ;H* k (0)) 

I emma B.2 : Let rand s, r 'and s 'be K -tuples of real numbers, such that for any k, 

1 <k<K,0<r k <r k and 0 s k < s k . We have for any V in <& h 

. K . ... s 4(r k " r k ) 2 (s i-’-s l ) 


(B. 12 ) H¥ II 3F'*\q) < c [L k=1 N(k) M(k) 


k -Y 


MV k (A k ;H}(0)) 


1/2 


Proof : The following inverse inequality, valid for s < s’ , is classical 
(B. 1 3) W€S m ( 0) , Hv|| s . e < C M s '- s ||v|| s0 ■ 

We derive immediatly (B.12) from (B.3) and (B. 1 3), 


Oorollaru B.2 : Let rand r'be two nonnegative real numbers, r < r'.We have for any ^ in 


P n , k (A) 0 S M (0) 


(B. 1 4) ||v||, „ < C ( N 2(r-r) * M r '- r ) [ I" s , ||v|| r % k ] 


1/2 


Proof : Once more, we deduce (B.14) from lemma B.2 and from the fact that H^(Q) coincides 
with H r (A; L 2 (0)) 0 L 2 (A; H^(©)). 


Now, for a given real number, we are interested in a-projection operator from H 0 ^(Q) 
onto^ h0 = fl H<5(A;L 2 (0)) . 

Theorem B.5 : There exists an operator TT,J from onto <& h0 verifying for any 

function V in H r ^(Q)H H<J(A ;L 2 (0) where r is a K -tuple of real numbers < 1 

(B. 1 5) || ^ — TT Jm/ 


h,n 


„ 2( 1 -r ) 2( 1 — r ) 

C[L k K -,[N(k) k + M(k) k ]|K|| 2 


1 1/2 


r k-°k 


Proof : The great difficulty in this Theorem relay on the fact that the degree in the Fourier 
direction are different in each Q k but the resulting approximation has to be globally in h’(O) 
which implies that at the boundary a k x©, the trace of the approximation must be in S y(k) with 
(j(k)s inf (M(k-l), M(k)). 
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Our proof will be decomposed into 3 steps, and we first assume that r k > 2. 

1) It is a standard result of the trace theory (see e.g. [12]) to note that if ^ belongs to 
H r '^(Q)fl Hq( A ;L 2 (0)) the first trace can be defined over each a k xO, more precisely, the 
mapping 


'*'-<* = \„e ■ ,<k<sl<+l > 

,1 


is continuous from H r '^(0)fl Hq(A ;l_ 2 (0)) fl L 2 (A ;H^(©)) into TT kz i H 


K+1 sup(r k _ r r k )-1/2 


(0). 


Besides, for any k, 2 < k K, there exists a continuous inverse mapping (see e.g. [12]) that 


associate to each 


sup(r, A ,r , ) - 1 / 2 

of H k -’ k (0), 


an 


element R k ( ij/ k ) of 


sup(r k-r r k ) 


#up(r k-iV 


H # R (A k-1 x0) x - (A k x0), whose first trace over a k x0 coincides with »j; k and 

whose first trace over a k _.,x0 and a k+1 x0 is zero. In particular the continuity of each R k that we 
can formulate as follows 

(B.16) V 4< k € H # (0) , ]] R k i|/ k ll 0k+ 1/2 < C II *k k ll 8k . 

implies that the mapping R that associate to each ^ of H r '^(Q)f] Hq(A;L 2 (0)) the element of 
L 2 (Q) defined as follows 


V k> 2 < k < K 1 , R^| Qk - [R k ^| akX 0]|Q k + [R k + i^|a k+1>< ©J|Q k 

R^i 


R Vr [R ^M ] i^ 


|0 


i w k 


K^|a K x0 ] |Q K 


satisfies R^ belongs to H s, *(Q)n Hq(A ;L 2 (0)) with s defined by 
(B.17) V k, 1 < k< K,$ k = inf( sup(r k _ 1 ,r k ), sup(r k ,r k+1 )) , 
where we have set r 0 = r K+1 = 0. Moreover, we also have 


(B.18) Vk, 


< k < K , (Hf - RV)| 0 € H,j(Q k ) 


2) It is an easy matter to find an element of <$& h0 that approximate ^ - Rt^. Indeed it suffices 
to take over each A k , (TT^ (k) ® n^ (k) )^ since this element vanishes over each elemental 
boundary. Moreover, it verifies 

(B.19) ||4/— (TT N(k) ® Tt M ( k) )^ Ili <fik <0( 11^- (TT N(k)® ^M(k))^ ^L 2 (A k ,H^(0)) 

+ ll^“ ( n N (k )< s>r[ M (k) )^ ; llHi(A .;l_ 2 (©)) - 1 ‘ 
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We write the first term 

ll^~ ^N(k)® ^M(k)^ ^L 2 (A k ;H^(0)) ^ II Od€> TT M )(M/— ( Tt N ® Id)^) |||_2^ A ; h^(0)) 

+ l?'-(id»nj)5'ii L 2 (A . H ^ (e)) . 

Now, we recall that TT^ commutes with the differentiation operators and satisfies for 0 < r < s 
[24]: 

(B.20) Vv<EH*(0), l|v-n;v|| r0 ^CM r - s ||v|| s0 . 

Using (B.5) for r = 0, (B.20) for r = s = 1 and (B.20) for r = 1 gives for s ^ 2, 


1 — r. 


ll i i'-(n' (u ®n' k) )i'|| L2( H , (e)) <c(N(k) k l ! H H v , (Ak . H ^ e)) 


k ’ V 


+M(k) k ll^ll|_ 2 ( A ; h\©))) • 


V ‘ V 


1 -r. 


1-r. 


Moreover , it is standard to note that, for r k 
(B.21) H^ k (Q k )cH rk_1 (A k ;Hi(0)). 

Hence, we obtain for s ^ 2 

(B.22) ||^-(n’ (k) ®n^ (k) )^|| L 2 (Ak . H i f(0)) <C(N(k) k +M(k) k ) ll^ll rk .n k 

In the same way, we write the second term 

II*- < n » l k)® n M<k)>*U h’CA.; L 2 (©)) < lk' d ® n m k ))(*-(T T » ( k)® | d)*) »H'(A. ;L 2 (©)) 


+ IIS'- (Iden^O^ II.. 1/ 




A .1 

i»i \ i\ / n \ / \ i ... 

.. V, . k , - 


Using (B.5) with r = 1 ,(B.20) with r □ s = 0 and (B.20) with r □ 1 gives 

»*- <n t!( k )® n M( kl >*H H ’< A k ; L 2 (©)) « C (N(k) '’ rk ll*llH rk (A k ; L 2 (©» 

*M(k)'‘ rk iifii H , (A H y' (0)) ) 


so that, by (B.21), 


1-r. 


1-r. 


(B.23) ll*-(n’ (k) ®n; (k) )*|| H i (A ;L 2 ( 0 )) <c(N(k) k + M(k) k mil rk .o k 


Finally, the inequalities (B. 1 9)(B.22)(B.23) imply 


1-r. 


1-r. 


( B .24) 11^- <TT,J (k) <8>njJ (k) )^ ll 1 . ftk ^C(N(k) k +M(k) k )|lH k .n k 

(cf, [14: Chap. I , Theoreme 5. 1 ]). 
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3) The same technique can be applied to the restriction of each x 0 ^k-i and • But 

here we know that the element R^u *0 *’ ias re 9 ulanty of ^ over Q k _, or Q k , We 

approximate R K MV by an element of IP n { k )( A k _ t ) x S m(k) over Q k _, and by an element 

of P n(k) ( A k ) x S m(k) over O k where we have set n(k) = inf (N(k- 1 ) , N ( k ) ) and 
m(k) = inf ( M ( k - 1 ) ,M ( k ) ) . Summing up the resulting approximations we derive an 
approximation n h ( v i') of Rst^ in <& h0 that satisfies 

1-sup(r r.) 1 -sup( r r.) 

(B.25) || Rt- - n h W)|U < C (n(k) +mCk) ) II RkVxS 'W,/ k ) 

r\ I 


1 — sup( r , r. ,) 1 -sup( r r ) 

+ n(k) k k+1 +m(k+l ) k k+1 


^ 1 a K+1 X © IUp(r k ,r k+1 ) 


Finally, we deduce that the polynomial Q h = ( TT ^ (k) <s> Tt^ (k) )^ + n h (\F) satisfies 
IK - Q h || 1(0 = II * + R¥ - Q h II,, Q < \\<m - (nj (k) ® n; {ti w ll,,, + ll ry - n h K) n,, n 
and the result (B. 15) follows, when all the p k are > 2, from (B. 1 7)(B. 1 8)(B.24) and (B.25) 
since 

|| V - TTK ||, n < K - Q h || 10 c [ L k=1 [ N(k ) 2<1 Pk> + M(k ) 20 rk> ] IK ll r 2 k ,n k ] 1/2 ■ 

The projection operator TT h is stable in the H^O) norm then we have 
K-TTh^|| 1fi <KI|,,, . 

The general result follows by using the main theorem of interpolation between Hilbert spaces of 

[H], 


We wish to obtain an error estimate for the interpolation operator. Let us denote 
by xr , f k = (i k ,j k ) , 0 < i k < N(k) , 0 < j k < 2M(k), the points (£ ik , 0 jk ) where 

0. k = -n + 2j k n/(2M(k)+ 1 ) . Let us defined the operator i h from C S (Q) into <5& h by 
(B.26) Vv€C’(0), Vk, 1 <k<K, VF k = (i k ,j k ) , I h v(x.-) = v(x.-) . 
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Theorem B.6 : Let rand s, r' and s' be K - tuples of real numbers greater than 1/2. We have 
for any function v in H r,s (Q) H ^ r,s (0) 

CB.27) llv- Ih v ll 0 ,(l <C [Lk = ,(N(k)' 2rk IMlV k (A t ;L 2 ( 0 )) +M(k) “ » v|I L 2 (A i; h\o)) 


V 
1-2r\ 


1-2s'. 


+ N(k) M(k) 




Proof : Let us denote by ijj the operator from C°(©) into $ M (0) defined by 
(B.28) V v € C*(0), V j , 0 j < 2M, I* v(ep □ v(0j) . 

Since over each Q k , i h is equal to i N(k) ® ijj[ (k) , we have 

ll v " I h v Ho,ft k < ( H v " (ld 0 I M(kP v ^L 2 (A k ; L 2 (0)) + " V " ^ X N(k) 0 ld)v HL 2 (A k ; L 2 (0)) 

+ || (Id - (Id ® I M(k)^^ v “ ^Nlk) 0 * d ^ v ^lL 2 (A k ; L 2 (0)) ' 
Then , we deduce the theorem from (B. 8) and the classical result (see [7; Theorems 1.2]), valid 
for s > 1/2, 

(B.29) V v € H ^( 0 ) , ||v- (Id® Im)v|| O0 ^CM _s ||v|| sQ . 

In the more simple case where N(k) and M(k) are independent of k, we have 
Corollaru B.3 : Let r be a real number > 1 .We have for any function v in H^(Q) and for 

< i o 

any \j > \ / 

(B.30) ||v- I h v|| 0ft ^C(M' r + (1+M _y N y )N 1/2 ' r ) ||v|| r0 . 

Proof : Using Theorem B.5 with r k = $ k = r k +$ k , independent of k together with (B.2 1 ) , yields 
l|v-I h v|l 0 ,o<C(M- r + N 1/2_r + M _r N 1/2-r+r ) ||v|| rfl . 

This inequality is valid for any r' = y > 1 /2. 

The end of this section being more related to our analysis than the previous general results, 
we shall suppose that N(k) and M(k) are independent of k; besides, it is clear, though tedious to 
write and read , how the general result would exist. 
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We want to estimate the difference. 

(B.3 1 ) E(f,w) = (f,w) - (f,w) hGL 

for any f in C°( 0) and w in P NK (A) ® S M . Let us denote by h* the couple ( N— 1 ,M) , we can prove 
Lemma B.3 : For any f in C°(Q) and anyw in P NK (A) <s> S M we have 
(B.32) | E(f,w) U C [ || f - TT h . f || 00 - || f - l h f || 0(0 ] || w || 0 0 . 

Proof : From (2. 1 5) 6L we deduce that 
E(f,w) = E(f-TT h , f , w) , 
hence, from (3.28) 

I E(f ,w) | < C [ || f - n h . f || 0 o + || l h f - n h . f || 00 ] || w || 0 
< C [ || f - TT f ||q q + || f - I h f Ho q ] || w || 0 f) . 


Using (B.9) and (B.27) yields immediatly 


Theorem B.7 : Let rand s, r’ and s' be real numbers greater than 1/2. For any w in 
P n k (A)<s>$ m and any f in H r (A;L 2 (©))n L 2 (A ; H s (0)) fl H r '(A ; H s '(0)) 

(B.33) | E(f ,w) | < C || w ||q (N llv llH r ( a • L 2 (0)) + M II v IIl 2 (a ■ H s (0)) 


N 


1/2 -r' M -s* 


H r ( A ; H s (0)) 


) 


Finally (B.9) and (B.30) gives 

Corollaru B.4 : Let r be a real number > 1. For any w in P NK (A) ® S M and any f in 
H r J,Ci) and for any y > 1/2 

(B.34) | E(f,w) | < C || w ||(M _r + ( 1 +M _t, N y )N 1/2 " r ) ||f || r<0 . 
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Figure caption 

Figure 1 . The error in the velocity u N = (u N ,v N ) and the pressure p N as a function of the total 
number of degrees of freedom ( Gauss Legendre Lobatto points ) in the x-direction when solving 
the discrete equation (4.12) corresponding to the test problem (4.13), The total interval 
A = ]-l ,1 [ is divided into two equal spectral elements of lengh A 1 = ]-1 ,0[ and A, = ]0,1 [, i.e. 
K = 2. Exponential convergence is obtained ( the plot in log-lin). 

Figure 2. The error in the pressure p N as a function of the total number of degrees of freedom 
( Gauss Legendre Lobatto points ) in the x-direction when solving the discrete equation (4. 1 2) 
corresponding to the test problem (4. 1 4). The error is given for x = 3 and for * = 5. The total 
interval A = ] — 1 , 1 [ is not divided into any subintervals, i.e. K = 1 . Algebraic convergence is 
achieved asymptotically ( the plot is log-log ) , although initially for small N faster convergence is 
achieved. 





□ o 
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